Time Optimal Unitary Operations 
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Extending our previous work on time optimal quantum state evolution [A. Carlini, A. Hosoya, T. 
Koike and Y. Okudaira, Phys. Rev. Lett. 96,060503 (2006)], we formulate a variational principle for 
finding the time optimal realization of a target unitary operation, when the available Hamiltonians 
are subject to certain constraints dictated either by experimental or by theoretical conditions. Since 
the time optimal unitary evolutions do not depend on the input quantum state this is of more direct 
relevance to quantum computation. We explicitly illustrate our method by considering the case of 
a two-qubit system self-interacting via an anisotropic Heisenberg Hamiltonian and by deriving the 
time optimal unitary evolution for three examples of target quantum gates, namely the swap of 
qubits, the quantum Fourier transform and the entangler gate. We also briefly discuss the case in 
which certain unitary operations take negligible time. 
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I. INTRODUCTION 



Time optimal quantum computation is attracting a 
growing attention [l], 0, 0] besides the more conven- 
tional concept of optimality in terms of gate complex- 
ity, i.e. the number of elementary gates used in a quan- 
tum circuit Q. The minimization of physical time to 
achieve a given unitary transformation is relevant for the 
design of fast elementary gates. It also provides a phys- 
ical ground to describe the complexity of quantum al- 
gorithms, whereas gate complexity should be regarded 
as a more abstract concept in which physics is implicit. 
Works relevant to the former subject can be found, e.g., 
in Q and 0, which discuss the time optimal generation 
of unitary operations for a small number of qubits us- 
ing a Cartan decomposition scheme and assuming that 
one-qubit operations can be performed arbitrarily fast. 
An adiabatic solution to the optimal control problem in 
holonomic quantum computation was given in [(|, while 
Schulte-Herbriiggen et al. 0] numerically obtained im- 
proved upper bounds on the time complexity of certain 
quantum gates. The present authors @] discussed the 
quantum brachistochrone for state evolution, i.e. the 
problem of finding the time optimal evolution and the 
optimal Hamiltonian of a quantum system for given ini- 
tial and final states. Nielsen ct al. Q proposed a cri- 
terion for optimal quantum computation in terms of a 
certain geometry in Hamiltonian space, and showed in 
@ that the quantum gate complexity is related to opti- 
mal control cost problems. Khaneja et al. suggested 
a geometrical method for the efficient synthesis of the 
controlled-NOT gate between two qubits with a special 
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Hamiltonian. 

In the standard quantum computation paradigm a 
whole algorithm may be reduced to a sequence of uni- 
tary transformations between intermediate states and a 
final measurement to read the result. In this paper we 
address the time optimality of each unitary transforma- 
tion, i.e., each subroutine. An example is the discrete 
Fourier transform in Shor's algorithm for factorization. 

In our previous work 0] , the quantum brachistochrone 
was formulated as an action principle for the quantum 
state in the complex projective space endowed with the 
Fubini-Study metric, and the Hamiltonian subject to cer- 
tain constraints. We obtained the time optimal state 
evolution and the optimal Hamiltonian by solving the 
Euler-Lagrangc equations. In the present work we ex- 
tend the methods used in Q and we describe the general 
framework for finding the time optimal realization of a 
given unitary operation. Roughly speaking, we replace 
the projective space representing quantum state vectors 
with the space of unitary operators. While the optimal- 
ity in the previous work depends on the initial state, it 
does not in the present case so that it is more directly 
relevant to subroutines in quantum computation, where 
the input may be unknown. This work should be useful 
not only for designing the efficient quantum algorithms 
and devices but also for deepening our insight into the 
true origin of the power of quantum computation. 

The paper is organized as follows. In Section [11] we 
introduce the problem by defining an action principle for 
the time optimal realization of unitary operations, un- 
der the condition of a Schrodinger evolution and of the 
existence of a set of constraints for the available Hamil- 
tonians, and we derive the fundamental equations of mo- 
tion. We discuss a typical class of the problem in Sec- 
tion lllll In Section llVl wc explicitly show how our formal- 
ism works via the example of a two-qubit system, which 
self-interacts by an anisotropic Heisenberg Hamiltonian 
depending on several control parameters. We derive the 
time optimal controls and the optimal time duration re- 
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quired to generate a SWAP gate, a 'qft' gate and an en- 
tangler gate. A system in which certain operations take 
negligible time is discussed briefly in Section [Vj Finally, 
Section IVII is devoted to the summary and discussion of 
our results. 



II. A VARIATIONAL PRINCIPLE 

Let us consider the problem of performing a given uni- 
tary operation or a quantum subroutine in the shortest 
time by controlling a certain physical system. Mathe- 
matically this is a time optimality problem of achieving 
a unitary operator Uf £ U (N) (modulo overall phases) 
by controlling the Hamiltonian H(t) and evolving a uni- 
tary operator U(i), where H(t) and U(t) obey to the 
Schrodingcr equation. Note that overall phases are physi- 
cally irrelevant for quantum evolutions. One immediately 
observes that there must be some constraints for H(t), 
because otherwise one would be able to realize Uf in an 
arbitrarily short time simply by rescaling the Hamilto- 
nian Q. Thus at least the 'magnitude' of the Hamil- 
tonian must be bounded. Physically this corresponds 
to the fact that one can afford only a finite energy in 
the experiment. Besides this normalization constraint, 
the available Hamiltonians may be subject also to other 
constraints, which can represent either experimental re- 
quirements (e.g., the specifications of the apparatus in 
use) or theoretical conditions (e.g., allowing no opera- 
tions involving three or more qubits). 

We then define the following action for the dynamical 
variables U(t) and H(i), 



S(U, H, A, Xj) := J dt [L T + L s + L c ] 



with 



Lj 1 

Ls : 
Lc 



L / <f.q-fr)(g)> 

V (HU,(l-Pu)(HU)) 



(1) 



(2) 

(3) 
(4) 



where we have introduced the Hilbert-Schmidt norm 
(A, B) := TrA^B and the projection Pu(A) := 
jfTr(AU^)U. The Hermitian operator A(t) and the 
scalars Xj(t) are Lagrange multipliers. The action term 
J Lxdt gives the time duration to be optimized and cor- 
responds to the action J[ds/v], where v is the velocity of 
the particle, in the classical brachistochrone. The metric 



dsj 



(dU,(l-Pu)(dU)) 



(5) 



is analogous to the Fubini-Study metric ds 2 FS = (di/)\(l — 
\t/))(ij}\)\dx(>) for the quantum state \ip) and is invariant 
under left and right global U(N) multiplications. 



The variation of Ls by A gives the Schrodingcr equa- 
tion 

i^L =H U, or U(t)=Te- i K Hdt , (6) 
dt 

where T is the time ordered product. This is similar 
to the case of the quantum brachistochrone for quantum 
states 0- On the other hand, the variation of Lc by Xj 
leads to the constraints for H, 



fj(H) = 0. 



(7) 



If we assume that the constraint functions fj (H) depend 

only on the traceless part of H, i.e. H := (1 — Pi)(H) = 
H - (TrH)l/N, thanks to the projection P v in © the 
action S is invariant under the U(l) gauge transformation 



U 



j U, Hi 



A, Xj » Xj, (8) 



where 9 is a real function. In the following we will con- 
sider the time optimal evolution of operators belonging 
to the group U(N)/ U(l) ~ SU(N). This is natural be- 
cause overall phases are irrelevant in quantum mechanics. 
To present our method in its simplest form, we have re- 
stricted ourselves to the case where the gauge degree of 
freedom is U(l). However, when there are quantum op- 
erations whose time duration is so short that it can be 
neglected, we will have a larger gauge group K. Such a 
case is discussed briefly in Section IVl 

We incidentally note here that, when the Hamilto- 
nian is time independent, the unitary operator actually 
evolves along a geodesic with respect to the metric ds^. 
This can be easily seen from ©, which implies 



£[(i-^)Of^)]=o, 



(9) 



the same equation as derived from the variation by U of 
the arclength J ds v . 

Let us now derive the other equations of motion. Be- 
fore taking variations of the action, it is convenient to 
rewrite Lx as 



Lt 



(Hil-P^H)) 



(10) 



where we have used the relation Pjj{A) = Pi(AU')U. 
Then the variation of S by H gives 



-L T .J 7 ZW1_ _ A + F = 0, 
(H,(1-P 1 )(H)) 

where we have introduced the operator 

dL c 



F := 



dH ' 



(11) 



(12) 



which plays an important role in the following. Using ^ , 
which implies Lt = 1, and recalling that (1 — Pi)(H) = 
H, one can rewrite (fTTj) as 



A = F 



H 



TriP 



(13) 
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Let us now take the variation of S by U. We first note 
that 

Tr As(^U^ = Tr D[A]U6rf (14) 

for any A up to a total time derivative, where D[A\ := 
^r u ^}- Tnc equation above holds because 8W - 



dt 

-WSUW and ^ 



W^W. Using (JTDJ) and (JUJl, 



one can easily calculate 5S/6U = to obtain 



D 



Lt 



(l_ Pl )(^tft) 



Pi)&)) 



iA 



0. (15) 



When the Schrodinger equation © and (Tl3"|) for A hold, 
one finds that the argument of D above is simply iF. 
We thus have D[F] = 0. Rewriting this, we obtain the 
quantum brachistochrone equation 



dF 
dt 



[H,F], 



This, together with the Schrodinger equation dp a nd the 
constraints ([7]), is our fundamental equation [111 ]. The 
quantum brachistochrone equation (|16[) seems universal, 
as it holds also in the case of time optimal evolution of 
pure and mixed fl2| quantum states [l3l |. In particu- 
lar, equation (JT5J) implies a simple conservation law, 



TrF' 



const., 



1,2, 



(17) 



In order to solve the quantum brachistochrone equa- 
tion ([To]) , one should first eliminate the gauge freedom 
([8]). The most natural gauge choice is to take H to be 
traceless, i.e., 



H = H. 



(18) 



This corresponds to choosing the unitary operator U to 
be an element of SU (N). Then, for a given operation Uf, 
the procedure to find the optimal Hamiltonian H and the 
optimal time duration T as follows: 

(i) specify the functions fj{H) which constrain the 
range of available Hamiltonians; 

(ii) write down the quantum brachistochrone equation 
(USD; 

(hi) solve ([To]) together with the constraints (J7J) to ob- 
tain H(t); 

(iv) integrate the Schrodinger equation ([6]) with 
U(0) = 1 to get U(t); 

(v) fix the constants in H (t) by imposing the condition 
that U(T) equals Uf modulo a global U(l), i.e., 



U(T) 



Uf, 



(19) 



where x is some real number. 

In essence, we have reduced the problem of finding the 
time optimal unitary evolution, for Hamiltonians sub- 
ject to certain constraints, to a set of first-order ordi- 
nary differential equations, which we call the quantum 
brachistochrone equation. Such an equation can always 
be solved in the general U (N) case, e.g. numerically. 



III. TYPICAL CLASS OF CONSTRAINTS 

Let us now discuss a typical and important class of 
constraints. We assume that the normalization condition 
for H, i.e., the finite energy condition, can be written in 
the form 



f(H) := \{TyH 2 -Nlo 2 )=Q, 



(20) 



where lo is a constant. Then the constraint part of the 
Lagrangian can be rewritten as 



L c = Xf(H) + L' c , 



(21) 



where A is a Lagrange multiplier and L' c is the sum of the 
other constraints. Therefore, from ([12]) . ([20[) and (|2"Tj) we 
obtain 



F = XH + F', 



(22) 



F(t) = U(t)F(Q)U\t). (16) where F' := 



dL' ( 
dH 



Multiplying by U from the right, 



using the quantum brachistochrone equation (|16p . the 
Schrodinger equation j6|) and (p~8|) . we have iXU + F'U — 
UF(Q). By formal integration, we get 



U = 



Texp 



F'dt 



exp (-£F(0)£| 



(23) 



The system becomes particularly simple if the con- 
straints for H are, except for the finite energy condition 
(|20[) . linear and homogeneous in H, namely, if 



L' c = Tr HF', 



(24) 



where F' = ^ Xjgj with gj G su(N), so that we have: 



TrgjH = 0. 



(25) 



Many problems in quantum computation or quantum 
control, including the example in the following section, 
fall into this subclass. Note that with the assumption 
F' does not depend on the Hamiltonian H explic- 
itly. 

We can easily show that A in ([22)1 is a constant. Choos- 
ing the gauge ([T5[) . we have 







dt dt 



(26) 



where the first equality follows from (TTu]) and the sec- 
ond one from the constraints ([20]) and ([25]) . Thus A is a 
constant, which can be chosen equal to one by a simple 
rescaling of F. From (|23|). we finally get 



U 



Texp N J F'dt 



exp(-iF(0)t), 



(27) 



while the Hamiltonian H(t) and the Lagrange multipliers 
Xj(t) are determined by ([To]) , i.e. 



dH 



dXj 
dt 



• XlF''' "-'A/- 



(28) 
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IV. EXAMPLE 

So far we have developed a general framework for find- 
ing the time optimal Hamiltonian. Let us now illus- 
trate our method by solving some specific examples ex- 
plicitly. What we consider is a physical system of two 
qubits represented by two spins interacting via control- 
lable, anisotropic couplings Jj(t) (j = x, y, z) and subject 
to local, controllable magnetic fields B a {t) (a = 1,2) re- 
stricted to the 2-direction. In other words, we choose as 
an example the following two-qubit Heisenberg Hamilto- 
nian, 



(29) 



Namely, the form 
guaranteed by 



29]) of the physical Hamiltonian is 



(32) 



]=x,y 



where Xjk (t) and (t) are Lagrange multipliers. 

Our task is to solve the quantum brachistochronc equa- 
tion (fT6|) , or ([28]) . Comparing the coefficients of the gen- 
erators of SU (4) on both sides, we find that the Lagrange 
multipliers \ xy and \ yx and the coupling J z are con- 
stants. Furthermore, the control variables B± and J± 
decouple from the others and we obtain, 



where cr j := a 3 ■ <g> 1, er| := 1 ® (Xj and <jj are the Pauli op- 
erators [14[ . In the standard computational basis labeled 
as |00), |01), 1 10), |11), the Hamiltonian ((29]) reads 





\-J z + B+ 








-J- 


H = 





J z +B- 


-J+ 








-J+ 


J z -B- 







-J_ 








-Jz~ 



B, 



(30) 



where we have introduced B±{t) := S 1 (i) ± B 2 (t) 
and J±(t) := J x (t) ± J y (t). By simply reordering 
the basis states as 1 00) , |11), |01), |10), the Hamiltonian 
can be rewritten as H 



Jz + B ± 
-J* 



-Jz — Ba 



= H + (B where H± := 

. We assume the finite energy 



condition (1201), i.e 



Bl + Bt + Ji + Jt 



2 J? = 2w 2 . 



(31) 



Then our problem is an example of the linear homo- 
geneous constraints discussed in the previous section. 



B±(t) = B 0± cos2( 7± i + ^±), (33) 
J±(t) = TS 0T sin2( 7T t + ^ T ), (34) 



X xy ± X yx are constants. 



where Bq± , tp± and 7± 

Let us now obtain U(t) by directly solving the 
Schrodinger equation ([6]) instead of using (|27p . Thanks 
to the block-diagonal form of the model Hamiltonian, 
the unitary evolution operator is also block diagonal, i.e. 
U = U+ ® U-. Using the Bakcr-Campbcll-Hausdorff for- 
mula (see, e.g. (l5j). in the permuted computational ba- 
sis the quantum evolution is described by the two decou- 
pled equations 



dV± 
i— — 
dt 



{TJz 1 ± l±v y + B ±a z )V±, 



(35) 



where V± := e^± t+ ^±>yU±. Solving ([35]) and going 
back to the original (non-permuted) computational basis, 
we finally obtain the optimal unitary evolution operator 
as 



e iJ * t (a 0+ + ia z+ ) 



U(t) = 



e Lht (a - + iu z -) 
s~ Uzt (—a y 



la-, 





e~ tJzt (a y - + ia x -) 

,-iJzt 



tot. 







'^{a y+ + ia x+ ) 



lJ **(a + - ia z+ ) 



(36) 



where we have chosen t/(0) = 1 and we have introduced 
the paramctrization 

oto±(t) ■= cos 7 ±i cosf2±£ + ff- smj±tsmQ,±t, 

a x± {t) := ±^sinQ ± tsin(7 ± t + 2V±) ! 

a y ±(t) := ±(smj±t cosO±t — cos7±tsinfi±i), 

a z ±(t) := -^sinO±tcos(7±t + 2V'±), (37) 

with ag ± + ol± = 1, and where fi± := J B%± + j±. 



The final step is to fix the coefficients J z , ao± (T) and 
a ±(T) (i- e - the constants -Bo±,7± an d ip±), the time du- 
ration T and the global phase x which realize the target 
([19]> . This is done imposing condition ([19]) with U(T) 
expressed via ([31)]) and ([3"7]) and Uf represented by the 
gate that we want to implement in a time optimal way. 
We now demonstrate this explicitly by a few simple but 
interesting examples. 

The swap gate: Let us assume that our target Ut is 
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the swap gate 



10 
10 
10 
1 



(38) 



which exchanges the states of qubits 1 and 2. Solv- 
ing (TT9")) by comparison of the matrix elements of ([3^)) 
and (f38|) and using (|37p . we obtain the following set 
of parameters: Bq + = 7 _ = 0, Bq-T = -|(1 + 2p), 
J*T = -f [1 - 2(p + g) - 4(?ti - n)], = f (1 + 2 9 ), 
and x = — j [1 — 2(p + q) + A{m + n)], where m,n,p 
and q are arbitrary integers and T is still to be de- 
termined. The zero values of B 0+ and 7 _, together 
with (f3"3")l and (fM| . imply that £?± and J± are con- 
stants and therefore, via (|3"D|) , that the optimal Hamil- 
tonian is time independent. The time optimal duration 
Tswap can then be found by imposing the constraint (|3"Tj) . 
which reads (4i£Zk*L) 2 = min TJl ,„, 1 ,^{2(l + 2p) 2 + [1 - 
2(p + <?) ^ 4(m — n)] 2 }. The solutions arc p = and 
q = —2(m — n) (or q = — 2(m — n) + 1), which lead to 
B X =B 2 = 0, J x = J y = 2J Z = (-1) 9+1 ^§, and finally 

give wT SWSP = ^7r and 



if 



V3 



" 1 








" 





-1 


2 








2 


-1 





. 








1 



(39) 



| (j> + q) + |(rn + n)], where again m, n,p, g, r and s arc 
arbitrary integers. As in the case of the U SWAP gate, the 
zero values of the parameters j± and eqs. p3|) and (|34|) 
imply that £?± and J± are constant and give a time in- 
dependent optimal Hamiltonian i7 SWAP , and consequently 
a geodesic evolution with respect to dsjj. Imposing the 

constraint (|3~lj) we obtain (^^) 2 = min m: „. p . ( j. I .{8(l + 
2p) 2 + [1 - 4(p + q) - 8(m - n)] 2 + 2(1 + 4r) 2 }', which is 
solved by p = r = and q = —2(m — n). This leads to 
B x = B 2 = J+ = J- = and finally gives the 

optimal time duration ll>T- QFT ' = -^-p-ft and the optimal 
Hamiltonian 



H 





~ 3 








0" 







-1 


4 





Vn 





4 


-1 







. 








-1 



(41) 



The entangler gate: As a last example, we want to find 
the optimal way to generate the entangler gate 



U rnt : — 



cos tp sin if 

10 

10 

— sin j) cos f 



(42) 



Since the Hamiltonian (|39[) is constant, due to © the 
time optimal evolution is along a geodesic for the metric 



dsfj. 



The 'qft' gate: Suppose now that we want to realize 
the slightly modified target operation U / given by 



10 
10 
10 
i 



This gate is important as it is essentially equivalent 
to performing a quantum Fourier transform (qft) over 
two qubits, i.e. U QFT = WiU' QFT 'Wi, where W\ is the 
Hadamard transform acting on qubit 1, i.e. W\ := 

" 1 and where the action of the qft on the 



l 

v^2 



1 



states of the two-qubit computational basis {\x) \ x = 
0,l,2,3} isgivenby@E/ QPT : \x) -> \ £!=q e^' 2 \y). 
The QFT is at the core of many quantum algorithms, such 
as the celebrated Shor's algorithm [l6[ for factoring in- 
tegers. If we can assume that the Hadamard transform 
takes negligible time, our methods generate the time op- 
timal Hamiltonian to obtain the target C/q F t- Following 
steps similar to those for the swap gate, we obtain the 
following set of parameters: Bq+T = (— l) r f (1 + 4r), 
B _T = |(1 + 2p), J Z T = -f [1 - 4(p + g) - 8(m - n)], 
7± = 0, 2>+ = vrr, 2^ = f (1 + 2q) and X = ~f [1 - 



where we choose the angle ip € [0,7r]. This gate, upon 
acting on the initial state |00), produces the (^-dependent 
entangled state cos^OO) — sin<^|ll). For example, when 
ip = 37r/4, 7r/4, this allows reaching the maximally en- 
tangled Bell states := (|00) ± \11))/V2. As usual, 
comparison of (|36|) and ([42]) leads to the following set of 
parameters: B 0+ T = (— l) r 7ry / j5 2 — (q — x) 2 , £?o- = 0, 



(40) J Z T = f (m - n), 7+ T = ir(q - x) and x = -f (m 



where x := f/n and, again, m,n,p,q and r are arbitrary 
integers. Imposing the constraint (|3ip we now obtain 
(^) 2 = min m ,„^. 9 {2[p 2 - (q - xf] + (m - ? i) 2 },which 
is solved by m = n and \p\ = —q = 1, and leads to 
B x = B 2 = (-l)V2wcos(2 7+ t + J x = Jy = 

(— l) p V / 2^sin(27 + t + We finally obtain the opti- 

mal, (/3-dependent wT ENT = tt^x(1 — x/2) and 



H(t) = ±V2u 



-cosp(t) smp,(t) 





s'mp(t) cos/i(i) 



(43) 



where n(t) := 2("/ + t+ip + ), 7 +(x) = uo(x—l)/yJx(\ — x/2) 
(see figure 1). In this case the Hamiltonian is time de- 
pendent and, therefore, the time optimal generation of 
the entangler gate does not occur along a geodesic. 
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FIG. 1: Time optimal realizations (|36p of the entangler gate 
(|42]) with if = k = 1, 2, 3, 4, from the identity (origin) in 
the a+-space. They are not geodesies. 



V. CASE WITH FAST OPERATIONS 

Let us now briefly discuss the quantum brachis- 
tochrone in the case where certain unitary operations 
have negligible time duration [19| . We assume that 
such operations (together with the unphysical global 
phase transformations) form a subgroup K of the group 
G = U (N) of all the unitary operations. We denote the 
Lie algebras of G and K by g and respectively. Note 
that the formulation in Section [IT] can be considered as 
the special case of K = U(l). 

In order to measure the time duration properly, 
we have to generalize the projection operator Pjj = 
P\{AW)U in ([2]) so that Pi is the orthogonal projection 
to 6 in g, namely, 

P 1 (A):=Y,(9j,A)g j , (44) 
j 

where {gj} is an orthonormal basis for 6. Apart from 
(HU), the action ((T|) and the Lagrangian terms ©-((1]) arc 
unchanged. Therefore, defining again H := (1 — Pi) (if), 
we can repeat the same argument and obtain the same re- 
maining equations of Section ILT1 and HLT1 In particular, we 
obtain the quantum brachistochrone equation (|16p and 
we can still follow the procedure (i)-(v) of Section HT1 

The terms Lt and L$ of the Lagrangian are now in- 
variant under the (in general non-Abelian) gauge trans- 
formation 

U^kU, H^kHk^+ikk\ 

K^kKk\ Xj^-Xj, (45) 

where k{t) £ K. Since F transforms as F i— > kFk^ un- 
der (|45|) . the quantum brachistochrone equation (fT6)) is 
always covariant, i.e., unchanged. Furthermore, if Lq 
is invariant under (|45p . the constraints are also gauge- 
invariant. 



One of the systems which is often discussed in quantum 
computation is that of n qubits in which the one-qubit 
operations take negligible time. This corresponds to the 
case G = U(2 n ) and K = U(l) ® SU(2)® n . Let S] 
be the subspace of g representing infinitesimal j-qubit 
operations. Namely, Qj consists of linear combinations of 
all the operators which are products of j Pauli operators 
and n — j identity operators: 

Qj := span t {<T°. ; ai < ■ ■ ■ < a 3 and lj = x, y, z}, 

J > 0, (46) 

where a rn represents the a m th qubit for m = 1, • • • , j and 
Oj^j" 3 is a generalization of er" appearing in (|29|1 . For 

example, cr^ = a x ® 1 ® a y ® 1 ® • • • £3> 1- Then we have 
6 = flo © 0i- Moreover, we can write ([44]) explicitly as 
Pi = £\ P« with 

pO)(A) : = x: E <C.; (*;'..::;" A ) 

(47) 

Note that each P^ is the orthogonal projection to Qj in 
0- 

Let us also assume that the infinitesimal operations 
including three or more qubits are not allowed in the 
Hamiltonian. This is the case of the linear homogeneous 
constraints discussed in Section (TTTJ with F' = Y^j=3^jj 
where 

/•; : E E m 

ai<---<a,j ii,-.. ,lj 

and A^ 1 ^ 3 are Lagrange multipliers. By choosing the 

gauge H = H we have Pj £ Qj , while from the constraints 

fj = with j > 3 we get P G 02- We find that the 
following commutation relations of the subspaces Qj of 
the algebra q hold: 

[&j,Qk\ = 8\j-k\+i ® 8\j-k\+s ® ■ ■ ■ ® Qj+k-i, j,k> 1, 
[0o,0j]=O, j>0, (49) 

where we understand that Qj — for j > n. 

In particular, the three-qubit case, n — 3, turns out 
to be simple and we can carry out the procedure in Sec- 
tion [IT] up to (iv) in general [19]. In fact, by (j4"9"|) we have 

[02,03] = 02, and since H £ q 2 and F' = P3 G 3 , the 
quantum brachistochrone equation (|16|) decouples into 
two equations 

P = -i[P,P'], F' = 0. (50) 

Thus we have F'(t) = P'(0) and H(t) = 
e lF ^* H(0)e~' lF so that we can drop the time or- 
dering in (|27p . and we finally obtain 

J7(t) = e ^'(0)t e -i(S(0)+F'(0))t_ (51) 
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A similar result was recently found in another setting 
[20j | . Although the three-qubit system is particularly sim- 
ple, the prescription (v) still remains technically involved. 
We postpone its full analysis to a future work. 

VI. SUMMARY AND DISCUSSION 

We have studied the problem of finding the time opti- 
mal evolution of a unitary operator in U (N) and the cor- 
responding time optimal Hamiltonian within the context 
of a variational principle. Our main result is an explicit 
prescription for finding the time optimal unitary opera- 
tion. Once the constraints for the available Hamiltonians 
are specified, the quantum brachistochrone equation can 
be immediately written down, and then the problem sim- 
ply reduces to obtaining its solutions. Our formulation 
is general, systematic, and does not rely upon any re- 
strictive assumptions, e.g., adiabaticity of quantum evo- 
lutions. We explicitly showed our methods and found 
the optimal Hamiltonian and the optimal time duration 
for three important examples of quantum gates acting 
on two qubits. The optimal Hamiltonians realizing the 
SWAP and QFT gates are time independent and, there- 
fore, the corresponding optimal unitary operators follow 
geodesic curves on the $J7(4) manifold endowed with the 
metric dsf,- This is not the case for the entangler gate 
(as expected for generic gates) where the optimal Hamil- 
tonian is time dependent and the time evolution of the 
corresponding unitary operator is not geodesic. We also 
discussed the quantum brachistochrone for unitary oper- 
ations in the case where there are operations whose time 
duration is negligible. 

This work is a natural extension of our previous anal- 
ysis of the time optimal evolution of quantum states in 
the projective space. The present formulation has direct 
relevance to quantum computation, since it gives the op- 
timal realization of subroutines for unknown input states, 
e.g. the discrete Fourier transform. On the other hand, 
the quantum brachistochrone for state evolution in Q 
may be viewed as a quantum computation for known ini- 
tial states, e.g. the transition from |00..0) to a certain 
entangled state [l6| in Shor's factorization algorithm. 

We should caution the reader that, in order to make the 
variational principle well defined, the action (p]) should 
be actually expressed as an integration over a parameter 
with fixed initial and final values. Since this does not 
affect our results, wc have omitted these details for sim- 
plicity. Furthermore, we note that, instead of j2]), any 
function of iUW and H which becomes constant upon 
using the Schrodinger equation would produce the same 
quantum brachistochrone equation (|16p . In this sense, 



the explicit expression of the metric in (|2|) does not af- 
fect our formulation. In a related work, the authors of 
[H , and [13] rephrased the problem of finding efficient 
quantum algorithms in terms of the shortest path in a 
curved geometry. Their goal was to obtain a bound on 
the number of gates required to synthesize a given tar- 
get unitary operation in terms of a cost function based 
on a certain metric in the space of Hamiltonians. By 
tailoring the form of such a metric they were able to ap- 
proximate the target unitary operation by a circuit of 
size polynomial in the distance from the identity. On 
the other hand, our point of view here is that the time 
complexity of an algorithm is of more physical relevance 
than its gate complexity (see also, e.g. Furthermore, 
although our result does not depend on the choice of the 
metric on U (N) , the bi-invariant metric ^ is the most 
natural. Also note that the simplest isotropic constraint 
(|20p does not provide any non-trivial bound to the gate 
complexity. In our framework the general relationship 
between time and gate complexity is still an open issue. 

Another point which we would like to emphasize is 
that, although what we treated here for the simplicity 
of exposition was the case in which the constraints are 
expressed as equality conditions for the functions fj (H) , 
there should be no conceptual difficulty in extending our 
variational methods to the more realistic case when sim- 
ilar constraints arc given in terms of inequalities (see. 
e.g., [13). 

Finally, we should note that the authors of [l[ , by using 
the Pontryagin maximal principle, also showed an opti- 
mal time dependent Hamiltonian as a particular solution 
to an equation which is similar to our quantum brachis- 
tochrone equation. In the two-qubit demonstration of 
our variational methods, we have obtained a general so- 
lution for the optimal Hamiltonian without attempting 
to match it to a prescribed NMR experiment, which was 
a main concern in Our formalism also naturally al- 
lows for the treatment of the more general and physical 
situation in which one-qubit local controls require a non- 
zero time cost. 
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